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I. INTRODUCTION 

Recent experimental discovery of Dirac and Weyl Q semimetals enhanced essentially the development of the 
interdisciplinary field of research related to the interaction between the condensed matter physics and the high energy 
physics. Together with ^He-A Q these materials are able to serve as an arena for the experimental investigation 
of various effects specific for the high energy physics because the low energy effective thec^ that describes 

fermionic quasiparticles in Dirac and Weyl semimetals has an emergent relativistic invariance [sl. [22 - 1^ . However, 
unlike the case of the superfluid Helium, the fermions in Dirac and Weyl semimetals are charged, and experience 
the external magnetic and electric fields. This facilitates, in particular, the investigation of effects related to chiral 
anomaly [^. [26l - [^ . 

Within the high energy physics there exist certain difficulties related to the inclusion into consideration of the 
gravitational background. In particular, there is an ambiguity in the expressions for chiral anomalyin the presence of 
gravity with torsion. Different methods of calculation give different expressions (see references , where those 

expressions are presented). At the same time, certain observable effects (mainly, in the condensed matter systems 
with emergent gravity) may be related to chiral anomaly in the presence of nontrivial geometry - for example, the 
appearance of Ko pni n force acting on vortices in superfluid helium Q , the appearance of the effects of anomaly in 
Weyl semimetals |T^ . ITTI . [TMl, [l8l4^ [ 2 ^ 113, (in particular, chiral magnetic effect 0) , etc. Therefore, the 
theoretical and experimental investigation of Dirac and Weyl semimetals is extremely promising since it should be 
able to resolve principal internal problems of the high energy physics. 

In the present paper we proceed our previous theoretical investigation of Dirac and Weyl semimetals [H, [l^, [13 ■ 
We extend the consideration of the 2D graphene (iol - IH . I46l - l48l| to the 3D Weyl semimetals. We consider the particular 
cubic tight - binding model [13, which is able to describe qualitatively a certain class of real materials. Although this 
model does not describe any particular semimetal, we do not exclude, that our results may also be relevant for the 
quantitative study of some real materials in a certain approaximation. The regular crystal, which is described by the 
given model represents a Weyl semimetal with two Fermi points. Each Fermi point hosts the two - component Weyl 
fermion. The chiralities of these two Weyl fermions are different. Such materials have more reach phenomenology, 
than the Dirac semimetals, where each Fermi point hosts the pair of Weyl fermions of opposite chirality (i.e. the 
massless Dirac spinor). 

We discuss the situation, when elastic deformations are present that result in the appearance of the emergent gauge 
field and emergent gravity. We calculate both emergent gauge field and emergent vierbein and express them through 
the tensor of elastic deformation. 

The paper is organized as follows. In Section |H] we describe the general algorithm for the calculation of the 
emergent fields in the presence of elastic deformations that was applied earlier to graphene. In Section IlIIl we recall 
the description of the unperturbed tight - binding model and demonstrate, that it indeed describes Weyl semimetal. 
In Sect 113 we consider the modification of hopping parameters of the tight - binding Hamiltonian resulted from the 
elastic deformations. In Sect. m we calculate emergent gauge field and emergent gravitations field in the presence of 
elastic deformations. In Sect. EH we describe the resulting low energy relativistic theory with the gauge field and the 
gravitational field. In Section rVHI we end with the conclusions. 


II. HOW DO THE EMERGENT GRAVITY AND EMERGENT GAUGE FIELD APPEAR IN WEYL 

SEMIMETALS 

Similar to the previous consideration of graphene [iol - ld^ we calculate the effective low energy action in the given 
tight - binding model of Weyl semimetal as follows: 

1. In the regular tight - binding model the positions of the Fermi points are calculated as the points in 

momentum space, where the one - particle Hamiltonian vanishes. The Hamiltonians are expanded around each 
Fermi point up to the terms linear in p — 

2. In the presence of elastic deformations the tight - binding model is modified. The terms in the Hamiltonian 

corresponding to the jumps between the adjacent lattice sites contain modified hopping parameters [13,113 ■ If 
those parameters would depend on the distance between the corresponding sites of the lattice onl y, t hen the 
relation between them and the tensor of elastic deformations would be similar to that of graphene [iol - liiL . 

Unli ke g raphene, however, the modification of hopping parameters in the 3D semimetals is more complicated 
(see ( 22 I I 43 I. We explore the modification of the hopping parameters following (22j| . 

3. It is assumed, that the variations of hopping parameters are weak. Therefore, we may consider space as a 
composition of regions such that within each region the hopping parameters do not depend on coordinates. 
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However, within each of such regions the hopping parameters are different for different directions of links. 

4. Within each area, where hopping parameters may be considered as constant we calculate the position of the true 

Fermi point K, where the one - particle Hamiltonian vanishes. This Fermi point is represented as K = + A, 

where A is interpreted as emergent gauge field. 

5. Next, the one - particle Hamiltonians are expanded around each true Fermi point K in powers of p — K. The 
coefficients of proportionality give the tensor of anysotropic Fermi velocity (within the given region) Q: 

H = flip, - Kf) (1) 

6 . The Hamiltonian in the whole sample is written as 

H = o-“/a(x) o (j), - - Ai(x)) = icr“ - A*(x)) + {p, - - A,(x))/* (x)) (2) 

Here the values of /^(x) and Ai(x) are calculated above within the regions being the small vicinities of space 
points X. The composition — Ai (x) is referred to as the floating Fermi point. 

7. Finally, the action for the fermions living near each Fermi point is written as [2^ 

S'=iy d‘^x\e\[^iel{x)<T^Vj'ii - [Dj^iel{x)a^'i>] (3) 

Here 

iVf,=iyf, + Af,{x) (4) 

is the covariant derivative (Aq = 0). It corresponds to the emergent 17(1) gauge field A^. By we denote the 

vierbein held, rjab is metric of Minkowski space. Internal 5'0(3,1) indices are denoted by Latin letters a, 6 , c,... 
while the space - time indices are denoted by Greek letters or Latin letters i,j,k,... The inverse vierbein is 
denoted by ef (it is assumed, that = <5]). The determinant of ef is denoted by |e|. By o-“ we denote Pauli 
matrices, and imply = 1 . 

The emergent vierbein is related to the tensor of anisotropic space dependent Fermi velocity as follows: 

/a = |e|e^, [eg]"^ = |e| = det^/^f, = 0, e° = 0, i, a = 1,2,3 (5) 

III. REGULAR TIGHT - BINDING MODEL 

Let us start from the tight - binding model of lattice Wilson fermions . It describes the 3D topological insulator 
with cubic crystal lattice and has the following Hamiltonian: 

Ho = i^V;(x + lj)(^it7V -r7°)^/>(x) + i^^(x)(m + 3r)7V(x) + (/i.c.) ( 6 ) 

x,i X 

Here the sum is over the positions x of the 3D cubic lattice and over j = 1,2,3. t and r are the hopping parameters 
while by 7 * the Dirac matrices (in chiral representation) are denoted. Vectors \j connect the nearest neighbor sites 
of the lattice. The given model describes Dirac fermion with mass m located in the vicinity of the the point p = 0 in 
momentum space. 

Extra term 


Hb = + (/i.c.) (7) 

X 

transformes the given model to the model that describes Weyl semimetal. In order to demonstrate this let us consider 
the Hamiltonian in momentum representation: 

H = ^ '0(k) ^ t sin(klj) — r cos(klj) 7 °^' 0 (k) + ^ '0(k)^(m + 3r) 7 ° + 63 7 ° 7 ^ 7 ®^' 0 (k) 

kj k 


( 8 ) 
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Here k is momentum of the quasiparticle. Next, let us define 

p. = (9) 

a 

(where a = jl^l is the lattice spacing), and 

F = m + r (1 —coskli) (10) 

i=l,2,3 

and apply the transformation ^( 7 *^— 7 *^ 7 ^)V”- This gives the transformed one - particle Hamiltonian (pj^ = (j>i,P 2 ))' 

jj ^ f vp±a + {F - b 3 )a^ -vps 

V -vp3 -vp±a - {F + 63 ) 0 -^ 

Here we take into account that the vectors \j form the triad of the cubic lattice. We introduce the dimensionless 
parameter 



V = ta, 


( 12 ) 


where a is the lattice spacing. It will be seen below, that this parameter has the meaning of the Fermi velocity in the 
xy plane. Dirac spinor has the form of the pair of two - component spinors 


ij} = 



(13) 


Let us denote energy by £(j)). Then we have: 


, _ vp3 , ^ 3 I 

^ £{p) + vp±a + {F + b 3 )a^ ^ 63 + ^ 


For small energies (close to the Fermi points) we arrive at 

■02 ~-- 


VP3 


r '01 


vp^a + {F + 63 ) 0-3 

and the Hamiltonian for the reduced two - component spinors receives the form 

= '!;pj_cr + a^{F - 63) + 


H, 


reduced 


vp±a + {F + 63 ) 0-3 


y2p^ + {F + 63)2 


2 2 
V P 3 


= , ?f+b3)0 + - ^ 3 ) + 


{F + 63) 22 


„2 I {F+b3y 

Px 4-— 


v'^p\ + {F + 63)2 

For . / -i <g; - it receives the especially simple form (up to the terms quadratic in p): 

y v^ + ^r “ 


'V Psi 


Hreduced — '^P_L<^ “h 


m + 63 - 


(m + 63)^TO + r ^ (1 - cosklj) - 63^ + 


2 2 
V P 3 


2 =1,2,3 


There are two Dirac points 


±K(°) = ( 0 , 0 , 


with 




(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


Let us denote 


^ = V 1 + ^ 3 ) = 


^ + ^irn + b3) 


( 20 ) 
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Notice, that in the original model of lattice Wilson fermions v = ta^ra^l, while m has the meaning of the physical 
mass. In this case the value of A would satisfy A ~ 1 in the limit a —>■ 0. In our case due to the smalness of u ~ 1/200 
the value of A may, in principle, differ from unity. However, we will see below, that rather natural assumptions about 
the components of the anisotropic Fermi velocity tensor will give A ~ 1. 

The Hamiltonians near these Fermi points have the form: 

H+ = vp±a + vX{p 3 - ~ ^ (21) 

V 03 + m 

and 


= 'op_Lcr - vX{p 3 + ^ = -a^(vp±a + vX{p 3 + (22) 

V 03 + m V V 03 + TO/ 

After the transformation ijj- —>■ of the two - component spinor at —we arrive at the model with the right 

- handed fermion at and the left - handed fermion at — 

The Fermi velocity is anysotropic and has the form 


f ( 0 ) = V 


(I 

0 


0 0 
1 0 
0 


(23) 


Notice, that the given low energy theory may be applied if |K^°)|a <C tt, so that the two Weyl points are close enough 
to each other and we may indeed use Eq. m instead of Eq. (HI). 

In the following we consider for simplicity the case, when 


ra'^ta = v<^ 1 . 


As it was mentioned, we assume, that the Fermi point value is much smaller, than tt/c, i.e. 


with 7 <C 1. Then 


^( 0 ) ^ 

^ V 


^7T fa 


fi = 2vX 


b3-m 

63 + TO 


= 2vX 


V 77r 


a (63 + to) 


In the typical Weyl semimetals the components of the anisotropic Fermi velocity tensor are of the same order. 
Therefore, 0(63 + to) ~ v^firX, which is consistent with the value of fermion mass to ~ v"fX^ (of the original 

topological insulator Eq. ([ 6 |)) much smaller than <C n/a. The value of A satisfies 

= \^irn + b3)\ < |^(m + 63 )| = |^(to + &3)| ~ < A (24) 

From here it follows, that A ~ 1. 


IV. MODIFICATION OF HOPPING PARAMETERS RESULTED FROM ELASTIC DEFORMATIONS 

In the presence of elastic deformations (5lj | the hopping parameters depend on direction and on the position in space. 
First let us consider the simplest model that relates hopping parameters with the tensor of elastic deformations. In 
this model the hopping parameter corresponding to the jump between the two sites x and x + Ij depends only on 
the real distance between these two sites given by r(x, Ij) = |u(x + Ij) — u(x)|, where u is the displacement vector. 
Therefore, we substitute the hopping parameter standing at the link (x,j) (summation over k and to is assumed): 

Pt IpJ^Ukm) = t{l - (25) 


and 


r -)■ r(l - Pr IpfUkm) = r{l - PrUjj) 


(26) 
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Here 

'^km — '^{dk'^m ~t“ O^Uk) 

is the linearized deformation tensor. 

Now let us consider the following complication. Let us assume, that the parameter t recieves extra correction due 
to the non - diagonal elements of the deformation tensor : 

ty yyi - Pt l^^lfukm) + tp^ ^ = ^7-^(1 - PtUjj) + tPf ^ (28) 

n^j n^j 


V. TIGHT - BINDING MODEL IN THE PRESENGE OF ELASTIC DEFORMATIONS 

A. Appearance of emergent gauge field 
The modified tight - binding model corresponds to the Hamiltonian 

+ ^ ■Uj„7°7” - r(l - 7 °) V'(x) 


H = 






+ i ^'0(x)^(m + 3r) 7 ° + 63 7 ° 7 ^ 7 ^^^(x) + {h.c.) 


Let us denote the trace of the deformation tensor by 

i=1.2,3 

Repeating the above steps we arrive at the one - particle Hamiltonian 

/ V Ei=1.2 - l^tUrt) + P't J2n^z 'U-rnPn) 

+ + r X)i=l, 2,3 1 - (1 - / 3 rM«)cOSkli -& 3 ^t 

-«Ei=1.2 CT*(k (1 - PtUzz) + P't UzuPn) 


(29) 


(30) 


H = 


-VP3{1- PtUzs) - VP[ J2n^3 '^inPn 


V 


-.*(1 - a« 3) - E„« »..A yy (1 _ + jya j 


(31) 

Here the product of momentum operator pk and the coordinate dependent function Uij{x) is defined as ^{pkUij{x) + 
Uij{x.)pk). Eq. (1511) has the form of Eq. (Hip with the substitution 

Pi^ Pz= Piil - PtUii) + / 3 ( ^ UinPn 

n^i 

m^F = m + r ^ j^l — (1 —/3rUii)cos kl^ 

i=1.2,3 

KiM + r ^ [^(1 -/3rUii)(l - coskli) 

i=l,2,3 

M = m + r PrU (32) 

The reduced Hamiltonian receives the form 

v3 


Hreduced = V^ PiCr'' + (M + 63 ) + T ^ (1 “ ;drMji)(l “ COSkl*) - 63 ^ + P, 


2 = 1,2 


2 d 2 

3 


2 =1,2,3 


E +^ 3 ) {m - 63 ) + V 


2 = 1,2 




(33) 


Here we neglected the terms proportional to ^r{M + 63 ) following the unperturbed model with A ~ 1. Recall that 
we consider the case, when 


ra<ta = r;<Cl 
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Besides, we assume, that the Fermi point value is much smaller, than tt/o, i.e. 


t-w _ 
iVg — 




= 77r/a, 7 <C 1 


The third component of the unperturbed Fermi velocity is given by 


= 2vX 


h-m „ , v"fTT 
= 2v A- 


63 + TO “ ’ a (63 + to) 

The requirement, that this component is of the same order as = v gives 

a (63 + to) ^ V -jttX 


Therefore 


2 2 

|^r(M + 63)1 = \ -^r{m + r/ 3 rU + b 3 )\ 


' 2 t 2 


m - 


■^ 3 ) +yul 


< 


'2t 


■h)] = 


'2v 


(to + 63)1 ^ 77rA <C A Ri 1 (34) 


This sequence of relations demonstrates, that we may indeed neglect ■§^r{M + & 3 ). The position of the floating Fermi 
point corresponds to the following values of P\ 


P = 


0 , 0 ,± 






The corresponding values of momentum are related to these values as follows: 


1 - /3iUii /3 (mi2 P'tUis \ [ Pi 

/3(uii 1 - f3tU22 P'tU23 P2 I = 

/ 3 ( u 31 f5'tU32 1 - l3tU33 J \ P3 


± 


0 

0 

—(m+r/ 3 ru )2 


The solution is (up to the terms linear in u): 


Pi \ 

/ 1 + PtUll 

-/3^U12 

-^'tUl3 

P2 

~ ± -P'tUll 

1 + hu22 

-l3'tU23 

P3 1 

\ -P'tU3i 

-P'tU32 

1 + l3tU33 


Vbj-r 


0 

0 

i" a 

_ mrprU \ 

h 'i — m? ' 


This gives the position of the floating Fermi point 

± K = ± ( - ,57,,. 1 + ft„,, - = ±(K<‘« + A). 

V bi — / V 


(35) 


(36) 


(37) 


(38) 


where the emergent gauge field has the following components 


All 

A 2 

^3 


-Pt'^23 




y/bf — m? 


(a 


'tU33 — 


mrPrU 's/bi — rn? 


bl- 




(39) 


2 A 

Notice, that ~ ~ ~ This is natural to suppose, that in the given semimetal all nonzero 

components of the emergent gauge field are of the same order. This requirement gives the natural constraint on the 
hopping parameter r: 


VJTT 

r ^ ~ TO 

a 


In this case both terms in the last row of Eq. (15^1) are of the same order. 



























B. Appearance of space depending anisotropic Fermi velocity 


The linearized Hamiltonians near the floating Fermi points have the form: 




H+=vJ2 + v(^P3- 


b3 + M 


(40) 


and 


= -a' 


Y V- ^ . Ib 3 -M\ 


(41) 


Therefore, tensor of anisotropic Fermi velocity receives the form 


/I 0 

0 1 


1 - PtUll /3(ui2 PtUi3 

P'fUii 1 - PtU22 Piu23 


^0 0 (1-|^)2^|^ j V Piu31 P[U32 l-PtU33 


that is 


f Ri r; 


( I - PtUll /3( mi 2 
I3[uii 1 - l3tU22 


P'tUi3 

P'tU23 


y P'tUsl PtU32 (1 - PtU33 - 

We may also rewrite this expression in the following way: 

/:«(^1 + - (Pt + Pt) J2 ^bSjbu,t), 

3 b 

where the unperturbed Fermi velocity tensor is given by Eq. 

The Hamiltonians at the (floating) Fermi points have the form: 

H+ = - K,) 

and 


= - 




(42) 


(43) 


(44) 


(45) 


(46) 


VI. RELATIVISTIC FERMIONS IN WEYL SEMIMETAL IN THE PRESENCE OF ELASTIC 

DEFORMATIONS 

Near each of the two Fermi points there is the left-handed/right - handed Weyl fermion. In the presence of elastic 
deformations the action for the right - handed Weyl fermion (incident at the Fermi point +K) has the form 24| : 


Sr = - d^x\e\['i>iel{x)a'’'Dj'i> - ['Dj'i>]iel{x)a^'i>] 


(47) 


Here 


iVf, = iVi^ +Af,{x) (48) 

is the covariant derivative. It corresponds to the 17(1) gauge field with the only nonzero component given by Eq. 
(15^ . The left - handed fermion belongs to the point — K and has the action 


Sl 


1 

2 


d‘^x\e\[^iei{x)a’’'Dj'^ — ['Dj'^]iel{x)a^'^] 


(49) 
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Here cr° = 1, cr® = —cr“ for a = 1, 2, 3 while 

= idfi, - A^^{x) 

The vierbein in the absence of elastic deformations is given by 

|err' = |e‘“>l=..r, = e<»'- = 0, 4»>» = 0, = 2'/»„ 

\ O3 “T TfL / 


(50) 


(51) 


where a, i,j,k= 1, 2, 3, while = vpfl has the meaning of the unperturbed anisotropic Fermi velocity with the the 
3x3 matrix 


/ = diag(z/ v = 2yj 


h-m 
63 + m 


(52) 


In the presence of elastic deformations the emergent vielbein is expressed up to the terms linear in displacement 
vector as follows: 


4 = /:(l + (A+/l.i;|^)i»)+/ff/i 


e^ = 0, eO = 0 




VF 

lei = uj’l i - 


-(1 + (/3t +/3r72-^)Fw) 
TF \ O 3 — 3 / 


(1 "^2)3“) 


a, J, j, k,n = 1, 2, 3 


(53) 


Here 


-/3tuii /3(ui2 /3(ui3 

U = ( -hu22 P'tU23 

Pt'^31 f^t'^32 { — PtU33 — 


(54) 


In another notations we have: 




(p' hrPrU 

SsSj ^2 _ ^2 




(55) 


The emergent gauge field is given by 


A, = 


/ N 


mrPrU N — mP- 


Ho = 0 


(56) 


(Recall, that u = '^u-) The analogy to graphene prompts that the values of /3t, /3(, and j3r similar to the Gruneisen 

parameter are of the order of unity. 

In the particular case, when Pt = —Pi^ Pr = 0 the expressions for the emergent fields are especially simple: 


Ai Pt^i3 

Uj = -PtUij 


VW^ 


V 


(57) 
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VII. CONCLUSIONS AND DISCUSSION 


In the present paper we expressed the emergent gauge field and emergent vierbein through the tensor of elastic 
deformations in the given model of Weyl semimetal. As follows from the general theory , the spin connection 

does not appear in the low energy effective field model of the fermionic quasiparticles. Therefore, we deal with a kind 
of emergent gravity with torsion, in which the emergent S'0(3,1) curvature vanishes. This is the case of the so - called 
teleparallel gravity (see [i^ and references therein). 

We follow the prescription for the modification of hopping parameters in the presence of elastic deformations 
proposed in Efij ] . The calculated expressions may be used for qualitative investigation of various effects in Weyl 
semimetals. In particular, the complete expression for the Stodolsky effect may be given. (This is the correction to 
Aharonov - Bohm effect due to torsion. In the topological contribution to this effect was given, while the total 
expression should contain the contributions proportional to the Gruneisen parameters /Sj,/?(,/3r-) Besides, having 
in hand the expression for the emergent vierbein we are able to calculate the Nieh - Yan term [23|, that, possibly 
enters expressions for the chiral anomaly in Weyl semimetals. Although we started from the model that does not 
refer directly to any particular crystal, we do not exclude, that our results may be relevant for the description of real 
materials in a certain approximation. This would open the possibility to check experimentally the appearance of Nieh 
- Yan term in the chiral anomaly, and, thus to resolve the ambiguity in the theoretical expression for chiral anomaly 
in the presence of gravity with torsion mentioned in the Introduction. 

The most promising development of the present study may be related to the similar consideration of the tight - 
binding models that describe real Weyl and Dirac semimetals (TaAs, Cd 3 As 2 , NaaBi). Those models should be more 
complicated because of the non - cubic crystal structure. Therefore, the direct extention of the method developed in the 
present paper to the consideration of these materials should be complicated technically (though obvious theoretically). 
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